Sudden death and long-living of entanglement in an ion-trap laser 
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The dynamical properties of quantum entanglement in a time-dependent three-level-trapped 
ion interacting with a laser field are studied in terms of the reduced-density linear entropy 
considering two specific initial states of the field. Allowing the instantaneous position of the 
center-of-mass motion of the ion to be explicitly time-dependent, it is shown that either sudden 
death of entanglement or survivability of quantum entanglement can be obtained with a specific 
choice of the initial state parameters. The difference in evolution picture corresponding to the 
multi-quanta processes is discussed. 

1 Introduction 

One of the most striking differences between classical and quantum correlations is the restricted 
capability of quantum states to share entanglement. The decay of entanglement cannot be re- 
stored by local operations and classical communications, which is one of the main obstacles to 
achieve the quantum computer [T|. Therefore it becomes an important subject to study the 
loss of entanglement El El El 13 CD] • Quite recently, by using vacuum noise two-qubit, entan- 
glement is terminated abruptly in a finite time have been performed [2J and the entanglement 
dynamics of a two two- level atoms model have been discussed jS] . They called the non-smooth 
finite-time decay entanglement sudden death. 

On the other hand, trapped atomic ions are an ideal system for exploring quantum infor- 
mation science. Recent advances in the dynamics of trapped ions (for a recent review, see e.g., 

have demonstrated that a macroscopic observer can effectively control dynamics as well as 
perform a complete measurement of states of microscopic quantum systems. With the reliance 
in the processing of quantum information on a cold trapped ion, a long-living entanglement in 
the ion-field interaction with pair cat states has been observed [5] . Also, experimental prepa- 
ration and measurement of the the motional state of a trapped ion, which has been initially 
laser cooled to the zero-point of motion, has been reported in [B]. 

In this paper we present an explicit connection between the initial state setting of the field 
and the dynamics of the entanglement. We give a condition for the existence of either entangle- 
ment sudden death or long-lived entanglement. In particular, a quantitative characterization 
of a general system of a three-level trapped ion interacting with a laser field is presented. We 
present various numerical examples in order to monitor the linear entropy and entanglement 



dynamics. The paper is organized as follows. In section 2, we consider a general class of a 
three-level system and obtain its solution. In section 3 we discuss the dynamics of the entan- 
glement with different initial states. Finally, we summarize the results and conclude in section 
4. 
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Figure 1: Energy-level diagram for a three-state A— type system interacting with a bimodal 
cavity field coupling the two ground states |1) and |3) to a common excited state |2) via a 
Raman transition. 

2 Model 

Now let us consider the Hamiltonian which describes a single trapped ion in a two-dimensional 
trap. Therefore, the physical system on which we focus is a three-level harmonically trapped ion 
with its center-of-mass motion quantized. We denote by ^ and ip\ the annihilation and creation 
operators and v 1(^2) is the vibrational frequency related to the center-of-mass harmonic motion 
along the direction x(y). The trapped ion Hamiltonian may be written as [TT| IT2l IT3] 

H = H + H int (t), 

2 3 

i=l i=l 

H int (t) = h%(x,t)S 12 + h^X 2 %(y,t)S 31 + hQl(x,t)S 21 + h%(y,t)S 13 . (1) 

We denote by Si m the atomic flip operator for the |m) — > \l) transition between the two 
electronic states, where Si m = \l){m\, (l,m = 1,2,3). 
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So far we have disregarded relaxations since we are interested in the dynamics for short 
times. Suppose the ion is irradiated by a laser field of the form 

Sfi(x,t) = 1 1 ® 1 exp[— i(kiX — flit)], 

%(y,t) = €2{lld2 h P2l2) eM-i(hy - mi (2) 

where e\ and €2 are the amplitudes of the two laser fields with frequencies fl\ and Q2 and 
polarization vectors pi and p2, respectively. The transition in the three- level ion is characterized 
by the dipole moment di and k{,i = 1,2 are the wave vectors of the two laser fields. We define 
the detuning between the atomic transitions and the fields as Ai = CJ21 — mifl\ and A2 = 

^23 - ^2^2- 

Therefore if we express the center of mass position in terms of the creation and annihilation 
operators of the two-dimensional trap namely 

x = Ax{^) x and y = Ay$ 2 + $ a ). (3) 

where Ax = (h/2vim) 1 ^ 2 = r\ x jh\ and Ay = (H/2v2m) 1 ' 2 = 772/^2 are the widths of the 
bi-dimensional potential ground states, in the x and y directions is called Lamb-Dicke pa- 
rameter describing the localization of the spatial extension of the center-of-mass in i— direction), 
and m is the mass of the ion. 

Making use of the special form of Baker-Hausdorff theorem the operator exp[ir]('ip 1 +'ip 1 )] 
may be written as a product of operators i.e. exp(i77(/0 +V>)) = exp fy] \ exp {ir^ J exp (irjip 



The physical processes implied by the various terms of the operator 

exp *, ft* ^ = exp =3- W E MjL. (4) 



. _ . n\ * — ' ml 

x ' n=0 m=0 

may be divided into three categories (i) the terms for n > m correspond to an increase in energy 
linked with the motional state of center of mass of the ion by (n — m) quanta, (ii) the terms with 
n < m represent destruction of (m — n) quanta of energy thus reducing the amount of energy 
linked with the center of mass motion and (iii) (n = m), represents the diagonal contributions. 
When we take Lamb-Dicke limit and apply the rotating wave approximation discarding the 
rapidly oscillating terms, the effective interaction Hamiltonian (0) takes the form 

Hint = h ll {t)£^\fy l )S 12 ^T 1 + h l2 (t)S^(^ 2 )S 2 ^T 2 + hrt(t)£M*$%)Sn^7 

W2(t)f*(tt a )wr, (5) 
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where 7j(t) is a new coupling parameter adjusted to be time dependent. The other contributions 
are rapidly oscillating with frequency v and have been disregarded. Note that in the Lamb- 
Dicke regime only processes with p — 0, 1 are considered, while in the general case, the nonlinear 
coupling function is derived by expanding the operator-valued mode function as 

Since S^fyi^j) depends only on the quantum number -0^-0^, in the basis of its eigenstates, 
j\n,i) = ni\rii), (n = 0,1,2,...), these operators are diagonal, with their diagonal elements 

(n|£^ ) (i/> i '0 i )|7i) is given by £^\rii) = -0.5ej(rij + fc)!) -1 ™*!!,*.^?) exp (-rfj2) where L k n {j]f) 

are the associated Laguerre polynomials. 

In what follows we obtain a general result regarding the solution to the time evolution 

operator. Now, we expand the time evolution operator in terms of the complete set of atomic 

operators as 

U(t) = exp ^-i J H int (r)d^j . (7) 

In the basis of the eigenstates, \ni,n2), we can find the elements of U(t) as &u(ni, ri2, t) = 
(i\U(ni,ri2,t)\j), where U(ni,ri2,t) = (ni, ri2\U(t)\ni, 712). Let us discuss the problem un- 
der consideration with Si = e, - = 0, and the initial conditions &u(ni, ri2, 0) = 1 and 
&ij(ni,ri2,0) — 0, (i ^ j). Under these conditions and after straightforward calculations, 
one can find an analytic time dependent solution in the following forms 



« ( +\ 1 / 7i(m + mi)! / /" \ 72(^2 + m 2 )\ 
&ii(n 1 ,n 2 ,t) = — I — cos U ni>n2 / 7 (r)drl+ — 

/V,n 2 L Ul - \ JO / n 2- 

6 22 (ni,n 2 ,t) = cos ^/i„ ijn2 ^ 7 (r)rfr^ , 

633K,^2,0 = -j S j cosU / 7(r)dr)+- 

/V,n 2 L n 2i V JO / 



ni! 



~ , ,x -i7i («i+ m i) ! • , 
6i 2 (ni,n 2 ,t) = 4/ j sinj/i 

A*ni,n 2 V 77,1 ' 



6i 3 (ni,W2,0 = 7^ (cos (V m , n2 / 7(r)^ - l) , 

A i ni,n 2 L V JO / J 



e.K.n.i) = -J2L j(J!l±^l sin L f\ (T)dT ). (8) 

The rest of 6jj(ni,n 2 ,t) = 712, t))* and the generalized Rabi frequency {i ni>n2 is given 
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by 



/7?£i 2 0i)(wi + mi)! 7l£f (n 2 )(n 2 + m 2 )! 
^ = V ^ + ^ • (9) 



We have thus completely determined the exact solution of a three-level trapped ion in the 
presence of time-dependent modulated function. Let's observe that putting in particular = 
1, 7(7-) = 1 and looking at the present solution (jSJ), we recover the time- independent three- level 
system discussed in [TTJ H2] . 

The (internal level) ionic dynamics depend on the distributions of initial excitations of both 
the field and the center-of-mass vibrational motion, given by p F (0) and p^(0), respectively. 
For instance, through a unitary evolution operator, the final state p(t) can be calculated in the 
following expression 

p(t) = U(t)(p A (0)®p F (0))U\t). (10) 

Having obtained the explicit form of the final state of the system p(t), we can discuss any 
statistical property of the system. We will turn our attention to the time evolution of linear 
entropy and entanglement, when the field is initially in the Fock state or coherent state. 



3 Entanglement dynamics 

The entanglement can be described by the linear entropy or the von- Neumann entropy [15J. 
The most prominent choice of pure state entanglement measures is the von-Neumann entropy 

[13 eg im nn 

s (pa(f)) = - tr {Pa[f) ^Pa(f)) , 

of the reduced density matrix, often simply called the entanglement E{jp) = S (pa_(f)) °f the 
pure state We work with the linear entropy which is convenient to calculate which 
is given by 

Sa (t) = l- tr A (p 2 A (t)) , (11) 

which ranges from for a pure state to 1 for a maximally entangled state and tr A denotes the 
trace over the subsystem A. The linear entropy is generally a simpler quantity to calculate 
than the von Neumann entropy as there is no need for diagonalization and can be considered 
as a very useful operational measure of the atomic state purity. 

In figure 2, numerical results for the time evolution of the linear entropy for an initial Fock 
state of the fields (|n, m), with n = m = 0) have been presented. In the typical experiments 
at NIST [20], a sm gle 9 Be + ion is stored in a RF Paul trap with a secular frequency along x 
of u/2tt ~ 11.2 MHz, providing a spread of the ground state wave function of Ax ~ 7 nm, 
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Figure 2: Example of the numerical calculation of Sa(0 as a function of the scaled time Xt, 
(A = 7j = 7 2 ). The parameters mj = 1, 7(r) = 1 and different initial states of the field, where, 
(a) Fock state with (n, = 0) and (b) coherent state with (rij = 0). 

with a Lamb-Dicke parameter of rj ~ 0.202. The two laser beams, with 0.5 W in each one, are 
approximately detuned A/2ti ~ 12 GHz, so that 7j/27r ~ 475 kHz. With these data we find 
€i ~ 0.01, so they can be considered as small parameters. The case of the effective vacuum is 
quite interesting where the linear entropy oscillates between zeros and a maximum value, in 
this case Sa (t) — 0.65 (see figure 2a). In fact the linear entropy attains the zero value (i.e., 
disentanglement) when the trapped ion is either in its upper or lower states (i.e., pure state) 
while strong entanglement occurs when the inversion is equal to zero. On other words, due to 
initial Fock state, the entanglement reaches its maximum value and drops to zero periodically, 
which opens the door for a possible application of the present model in constructing a quantum 
logic gate. 

The immediate question now is, if different initial states of the field are considered, is the 
periodic behavior of the linear entropy and zeros entanglement for such states still exist? To 
answer this question we make use of a coherent state as an initial state of the field and find a 
general entanglement feature, captured in equation (JTTj) and illustrated in figure 2b. Once the 
initial state of the field is considered to be a coherent state the situation is changed drastically 
(see figure 2b). It is obvious that the time evolution of linear entropy behaves as that of standard 
single-photon Jaynes-Cummings model and oscillates irregularly with the time. At the early 
times the linear entropy from zero evolves to its local maximum value ( ~ 0.96). In this process, 
the three level-trapped ion and the fields are always entangled. Although increasing the mean 
photon number leads to strong entanglement (maximum value of entanglement), however the 
maximum value of the entanglement also varies and occurs for some short period of time. This 
indicates that in a regime where coherent state is considered, the underlying states are highly 
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entangled. 




Figure 3: The same as figure 2a but in this figure we consider the modulated function to be 
time-dependent of the form (JUJ), the initial time is to = — lOr and different values of rii, where 
(a) rii = and (b) rii = 15. 

Let us consider the modulated function j(t) to be time-dependent of the form [2*T| 12*2] 

7(t) = sech 

In this form the coupling increases from a very small value at large negative times to a peak 
at time t = 0, to decrease exponentially at large times. Thus, depending on the value of r and 
the initial time to, various limits such as adiabatically or rapidly increasing (for to < Xt < 0) or 
decreasing (for < At < to) coupling can be conveniently studied. This allows us to investigate, 
analytically, the effect of transients in various different limits of the effect of switching the 
interaction on and off in the ion-field system. The vanishing of the interaction at large positive 
times leads to the levelling out of the inversion. It should be noted that the time dependence 
specified in (J12)) is one of a class of generalized interactions that offered analytical solutions. 

In figure 3 the changes in linear entropy vs the dimensionless quantity Xt is plotted when 
the modulated function is taken to be time-dependent as in equation (JT2J). An intriguing result 
found in figure 3a, where linear entropy is plotted with an initially Fock state of the fields 
(\ni,n 2 ), with ri\ = n 2 = 0) with a time-dependent modulated function, showing clearly the 
sudden death of entanglement at Xt ~ 10.2. A remarkable property of such initial state setting 
is that entanglement can fall abruptly to zero for a very long time and the entanglement will 
not be recovered i.e. the state will stay in the disentanglement separable state. On the other 
hand, we notice that the long-living entanglement can be obtained with large values of the 
initial Fock state numbers, such as \n\ = 15, n 2 = 15), (see figure 3b). Therefore, the initial 
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Fock states placed at this point is a suitable choice to investigate entanglement dynamics for 
different initial number of photons for the fields. It's not surprising to find that the number 
of oscillations is increased for higher mj. At the period —10 < Xt < 10, the slight difference 
lies on the number of oscillations only, while for the later times (say At > 10), the situation 
becomes completely different. 
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Figure 4: The same as figure 3 but in this figure we consider an initially coherent state field 
with an average photon number of Hi = 10 and various value of the number of quanta, where 
(a) rrii = 1 and (b) m; = 2. 

The above results and connections are very intriguing, and lead us to ask what is the role 
played by the initial state in obtaining these associated phenomena of the entanglement. In 
order to answer to this question, we consider different initial state in figure 4. This figure 
shows the linear entropy with an initially coherent state field with an average photon number 
of Hi = 10 for various number of quanta. The initial time is t = — lOr, so the interaction starts 
at a fairly low value, peaks and then drops off again. It is interesting to mention here that, as 
time goes on long-livid entanglement is observed. As a particular but striking enough example 
we have considered the same value of the system parameters which have been considered in 
the time-independent case. We have analyzed the long-lived entanglement by considering a 
multi-photon interaction (mi = m<i = 2) in figure 4b. This case is similar to a situation where 
both rrii — m 2 — 1, because both linear entropies in figure 4a and 4b rise and lower together 
although they are not equal. The only difference between the two cases is that the maximum 
long-living entanglement for one photon case is S(t) ~ 0.69, but in the two-photon case is 
S(t) ~ 0.62. Also, figure 4b demonstrates that the linear entropy peaks show a lowering of the 
local maximum at the interaction period — 10 < Xt < 0. 

All these results confirm the possibility of a practical observation of time- dependence of 
the modulated function effects for creating sudden death or long-lived entanglement. Based on 
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such sensitivity and some other evidence, we suspect that the analytical results presented here, 
could be attained for different configurations of the three-level systems. 

4 Conclusions 

Summarizing, we have investigated the dynamics of quantum entanglement for a trapped ion- 
laser field interaction. An explicit expression is given for a time-independent case and compared 
with previous studies. Through a three-level trapped ion system we have shown that the 
commonly assumed initial state setting may affect entanglement in a very different manner. 
This study reveals that the time-dependent modulated function can be used for generating either 
entanglement sudden death or long-lived entanglement depending on a proper manipulation of 
the initial state setting. We hope the presented results can be useful for the ongoing theoretical 
and experimental efforts in multi-levels particles interaction. Hence, despite the considerable 
progresses on which we have reported here, a panoply of challenging open questions awaits 
solution, what simply reflects the decoherence effect as well as the cavity decay or atomic 
decay. 
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